Effect of neutron composition and excitation energy on the primary fragment yield 

distribution in multifragmentation reaction 
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The isotopic properties of the primary and secondary fragment yield distribution in the multi- 
fragmentation of 58 Fe + 58 Ni and 58 Fe + 58 Fe reactions are studied with respect to the 58 Ni + 
58 Ni reaction at 30, 40 and 47 MeV/nucleon. The reduced neutron and proton densities from the 
observed fragment yield distribution show primary fragment yield distribution to undergo strongly 
secondary de-excitations. The effect is small at the lowest excitation energy and smallest neutron- 
to-proton ratio and becomes large at higher excitation energies and higher neutron-to-proton ratio. 
The symmetry energy of the primary fragments deduced from the reduced neutron density is signifi- 
cantly lower than that for the normal nuclei at saturation density, indicating that the fragments are 
highly excited and formed at a reduced density. Furthermore, the symmetry energy is also observed 
to decrease slowly with increasing excitation energy. The observed effect is explained using the 
statistical multifragmentation model. 

PACS numbers: 25.70.Mn; 25.70.Pq; 26.50.+X 



I. INTRODUCTION 

In multifragmentation reactions an excited nuclei ex- 
pand into vacuum and decay into various light and heavy 
fragments 0, 0, 0, H ■ These fragments when emitted are 
highly excited and neutron-rich, and subsequently un- 
dergo de-excitation to cold and stable isotopes. Similar 
hot nuclei are also produced in the interior of a collaps- 
ing star and subsequent supernova explosion The 
production of these nuclei depends on their internal exci- 
tation and is sensitive to the symmetry energy part of the 
binding energy jy] . A slight decrease in the symmetry en- 
ergy coefficient and subsequent de-excitation can signif- 
icantly alter the elemental abundance and the synthesis 
of heavy elements Q. In multifragmentation reactions 
the measurement of fragment isotopic yield distribution 
can provide important insight into the symmetry energy 
and the decay characteristics of these nuclei. 

Experimentally, the determination of fragment isotopic 
yield distribution is often influenced by the complex de- 
excitation of the primary fragments to their finally ob- 
served secondary yield distribution. Theoretical calcula- 
tions of the yields after secondary de-excitations require 
an accurate accounting of the feeding from the particle 
unstable state and are subject to uncertainties in the lev- 
els that can be excited, and the structure effects that 
govern their decay. 

In this work, we study the primary fragment yield dis- 
tribution and the associated symmetry energy through 
relative reduced neutron and proton densities at vari- 
ous excitation energy and isospin (N/Z) of the multifrag- 
menting source. It has been shown pj that the relative 
reduced neutron and proton densities can provide an al- 
ternate way of studying the isotopic yield distribution, 
where some of the uncertainties mentioned above can be 
avoided by taking the ratios of the yield in two different 
reactions that differ only in their isospin (N/Z) content. 



Also, the reduced neutron density is directly related to 
the symmetry energy through a scaling relation as dis- 
cussed in 0, 0111 EI- 

In this work we show that the primary fragment 
yield distribution is strongly affected by the secondary 
de-excitations at higher excitation energies and isospin 
(N/Z) of the fragmenting system. The symmetry en- 
ergy of the primary fragments deduced from the scaling 
parameter is considerably lower than that for the satu- 
ration value and decreases with increasing excitation en- 
ergy. The observed effect can be well explained by the 
statistical multifragmentation model of Botvina et al. [|| . 



II. EXPERIMENT 

The measurements were carried out at the Cyclotron 
Institute of Texas A&M University (TAMU) using beams 
from the K500 Superconducting Cyclotron. Isotopi- 
cally pure beams of 58 Ni and 58 Fe at 30, 40 and 47 
MeV/nucleon were bombarded on self-supporting 58 Ni 
(1.75 mg/cm 2 ) and 58 Fe (2.3 mg/cm 2 ) targets. Six dis- 
crete particle telescopes placed inside a scattering cham- 
ber at laboratory angles of 10°, 44°, 72°, 100°, 128° and 
148° were used to measure the fragments from the reac- 
tions. Each telescope consisted of a gas ionization cham- 
ber (IC) followed by a pair of silicon detectors (Si-Si) and 
a Csl scintillator detector, providing three distinct detec- 
tor pairs (IC-Si, Si-Si, and Si-Csl) for fragment identifica- 
tion. The ionization chamber was of axial field design and 
was operated with CF4 gas at a pressure of 50 Torr. The 
gaseous medium was 6 cm thick with a typical threshold 
of ~ 0.5 MeV/A for intermediate mass fragments. The 
silicon detectors had an active area of 5 cm x 5 cm and 
were each subdivided into four quadrants. The first and 
second silicon detectors in the stack were 0.14 mm and 
1mm thick, respectively. The dynamic energy range of 
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the silicon pair was ^16-50 MeV for 4 He and ~ 90 - 
270 MeV for 12 C. The Csl scintillator crystals that fol- 
lowed the silicon detector pair were 2.54 cm in thickness 
and were read out by photodiodes. Good Z identification 
was achieved for fragments that punched through the IC 
detector and stopped in the first silicon detector. Frag- 
ments measured in the Si-Si detector pair also had good 
isotopic separation. Further details can be found in ref. 

mm 

The present study is carried out for fragments detected 
at 44 degrees in the laboratory system, which corresponds 
to the center of mass angle « 90 degrees. The fragments 
detected at this angle originate predominantly from cen- 
tral events which was further assured by gating on the 
measured neutron multiplicities. A possible admixture 
of the projectile/target fragmentation processes will nev- 
ertheless be accounted for in the following analysis. Fig- 
ure 1 shows the experimentally determined yield distri- 
butions for the lithium (left panels) and carbon (right 
panels) isotopes measured in 58 Fe + 58 Fe (triangle sym- 
bols), 58 Fc + 58 Ni (circle symbols), 58 Ni + 58 Fe (square 
symbols) and 58 Ni + 58 Ni (star symbols) reactions at the 
three different beam energies. The data show the high- 
est yield for the neutron rich isotopes in the 58 Fe + 58 Fe 
reaction (triangles), which has the largest neutron-to- 
proton ratio (N/Z), compared to the 58 Ni + 58 Ni reaction 
(stars), which has the smallest neutron-to-proton ratio. 
The yields for the reactions with the intermediate value 
of the neutron-to-proton ratio, 58 Ni + 58 Fe (squares) and 
58 Fe + 58 Ni (circles), lie in between. One also observes 
that the difference in the yield distribution for the four 
reactions decreases with increasing beam energy. As will 
become evident in the following sections, this is a con- 
sequence of the secondary de-excitation of the primary 
fragments which becomes dominant for systems with in- 
creasing neutron-to-proton ratio and excitation energy. 



III. THE REDUCED NUCLEON DENSITIES 

Using the experimentally observed yield distribution, 
the reduced neutron and proton densities were extracted 
as follows. In a Grand-Canonical approach for the mul- 
tifragmentation process (see e.g. 

m El El El), the 
fragment yield for an isotope with neutron number N, 
and proton number Z (mass number A = N + Z), is 
given as 

Y(N, Z) oc Vp»p z p Z N<z {T)A^e B ^ z V T (1) 

where V is the volume of the system and p n oc e M "/ T 
and p p oc e^p/ T , are the primary 'free' neutron and proton 
densities, respectively. The exponents p n and p p are the 
neutron and the proton chemical potentials, and Z^,z{T) 
is the intrinsic partition function of the excited fragment. 
The quantity B(N, Z), is the ground state binding energy 
of the fragment and T, is the temperature. In the above 
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FIG. 1: Relative yield distribution of the fragments for the 
lithium (left) and carbon (right) isotopes in 5S Fe + 58 Fe, 58 Fe 
+ 58 Ni, 58 Ni + 58 Fe and 58 Ni + 58 Ni reactions at three dif- 
ferent beam energies. The different symbols are for the four 
reactions as discussed in the text. 

formula, the effect of the Coulomb interaction on frag- 
ment yield is often neglected. By introducing p n and p p , 
the actual isotope yields are reduced to an approxima- 
tion appropriate for the thermodynamical limit at high 
excitation energies E3- Hence, they are referred to as 
the 'reduced' densities. 

As discussed in the introduction, taking the ratios of 
the fragment yield from two different systems which dif- 
fer only in their isospin (N/Z) content reduces the uncer- 
tainties in the quantities shown in Eq. I. For the present 
analysis all ratio's were taken with respect to the least 
neutron-rich 58 Ni + 58 Ni reaction. The isotopic yield 
distribution of the fragments in terms of relative reduced 
neutron density is then given as, 



Y(N + k, Z)/Y Nl+m {N + k, Z) 
Y(N, Z)/Y m + Ni (N, Z) 




where fc, correspond to the different isotopes used to 
determine the double ratio, and Y Nl+Nl is the yield for 
the 58 Ni + 58 Ni reaction. A similar expression for the 
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FIG. 2: Experimentally determined relative reduced neutron 
and proton densities as a function of excitation energy for the 
58 Fe + 58 Ni and 58 Fe + 58 Fe reactions with respect to the 58 Ni 
+ 58 Ni reaction. The different regions bounded by the solid 
and the dotted lines are the statistical model calculations as 
explained in the text. 



relative reduced proton density from the isotonic yield 
ratios can also be written as 



Y(N, Z + k)/Y m+m (N, Z + k) 
Y(N, Z)/Y m + Ni (N, Z) 



m- 
Pp 



■Ni 



(3) 



Figure 2 (symbols) shows the relative reduced neutron 
and proton densities as a function of excitation energy 
(which correspond to different beam energies as discussed 
in the next section) of the fragmenting source for the 58 Fe 
+ 58 Ni and 58 Fe + 58 Fe reactions. The various symbols 
in the figure correspond to the densities obtained using 
k = 1, k = 2 

and k = 3, for various isotopes and isotones. The figure 
shows a steady decrease in the reduced neutron density 
and an increase in the proton density with increasing ex- 
citation energy. The effect is stronger for the 58 Fe + 58 Fe 
reaction, which has the highest isospin (N/Z) content. 



IV. COMPARISON OF THE REDUCED 

NUCLEON DENSITIES WITH THE 
STATISTICAL MULTIFRAGMENTATION 
MODEL CALCULATION 



In order to study the effect of secondary de-excitation 
at various excitation energies and isospin content (N/Z) 
of the fragmenting system, a comparison of the observed 
reduced neutron and proton densities with the statisti- 
cal multifragmentation model (SMM) of Botvina et at, 
jf| was carried out. The initial parameters for the SMM 
calculation such as, the mass (A), charge (Z) and exci- 
tation energy (E* ) of the thermal source were estimated 
by simulating the initial stage of the collision dynam- 
ics using the BNV model calculation [l7|. The param- 
eters obtained for the 58 Fe + 58 Ni reaction are, A s w 
111, Z s k 52 and E* w 5, 7 and 9.4 MeV/nucleon for 
the three beam energies studied. The N/Z ratios of the 
sources in all cases were not very different from those 
of the interacting nuclei. These results were obtained 
at a time around 40-50 fm/c after the projectile fuses 
with the target nuclei and the quadrupole moment of the 
nucleon coordinates (used for identification of the defor- 
mation of the system) nearly goes to zero. Also, the form 
of single source in the calculation oscillates throughout 
the dynamical evolution and does not undergo any dy- 
namical disintegration, assuring that the dynamical flow 
produced in these reactions are insignificant for the en- 
ergies studied. 

The source parameters were also compared with the 
experimental data of Hudan et al. |l8j where a nearly 
symmetric central collisions of Xe + Sn was studied 
and an excitation energy of w 5 MeV/nucleon at 32 
MeV/nucleon, and « 7 MeV/nucleon at 50 MeV/nucleon 
was obtained. It was shown in this work that the SMM 
calculation describes the disassembly of these systems 
and the excitation energies of the observed primary frag- 
ments quite well. A radial flow of « 2 MeV/nucleon was 
however identified in this case for the 50 MeV/nucleon 
beam energy. In the present case, the systems are twice 
as small and stable. Further, it is obvious from various 
measurements (see e.g. that smaller systems accu- 

mulate smaller radial flow than larger systems. Conse- 
quently, a large part of the energy becomes available in 
the form of thermal excitation. 

To account for the possible uncertainties in the source 
parameters due to loss of nucleons during pre-equilibrium 
emission and the sensitivity of the final results toward 
the source size, the calculations were also performed for 
smaller source sizes of A s = 58 with Z s = 26, 27, and 
28. The freeze-out density in the calculation was as- 
sumed to be 1/3 of the normal nuclear density and the 
de-excitation of the hot primary fragments was carried 
out via evaporation, fission or Fermi breakup 0] 

Figure 3 shows the mean characteristics of the hot 
primary fragments in the freeze-out volume calculated 
from the statistical multifragmentation model. Shown 
are, their relative yields, internal excitation energies and 
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FIG. 3: Characteristics of the hot primary fragments ver- 
sus their mass number A, produced in the freeze-out volume 
during multifragmentation of the thermal sources, with mass 
number A s — 116, and charges Z s — 52 and Z s — 56, at exci- 
tation energies of E* — 6 and 10 MeV/nucleon (see notations 
in the figure). Panels from top to bottom: neutron to proton 
(N/Z) ratio; difference of N/Z ratios between neutron rich 
(Z s — 52) and neutron poor (Z s — 56) sources at different 
E*; internal excitation energies; relative mass yields. 

N/Z ratios as a function of their mass. At the lowest 
excitation energy E* = 6 AMeV, the primary fragments 
are characterized by an exponential-like mass distribu- 
tion (bottom panel). With increasing excitation of the 
source the internal excitation energy of the fragments 
increases and the source breaks up into smaller frag- 
ments. For sources with different neutron-to-proton ra- 
tios (N/Z), the primary fragments are also produced with 
varying N/Z ratios depending on the source excitation 
energy (top panel). 

The characteristics of the hot primary fragments 
however, change significantly after the secondary de- 
excitation as shown in Figure 4. Both the mass and the 
neutron content of the cold secondary fragments decrease 
considerably lowering the N/Z ratio 

due to a higher evaporation probability of neutrons 
relative to protons. The N/Z of the light fragments vary 
with Z reflecting their shell structure. However, the qual- 



FIG. 4: Same as in Fig. 3, but for the cold fragments produced 
after the secondary de-excitation of the primary fragments. 



itative difference between the N/Z ratios of the hot frag- 
ments (middle panel) produced by the sources with dif- 
ferent isospin remains the same even after the secondary 
de-excitation, though it becomes smaller. Moreover, one 
can see the solid and the dotted lines in the middle panel 
to be more distinguished than in the Figure 3. The aver- 
age difference between the cases with the two excitations 
increases from 0.0051 to 0.0168 for the fragments A = 6 
- 20. This means that the difference in the N/Z of the 
fragments produced between the neutron rich and the 
neutron poor sources decreases more rapidly for high ex- 
citation energies as a result of more intensive secondary 
de-excitation. 

Figure 2 presents the comparison between the experi- 
mentally (symbols) determined relative reduced neutron 
and proton densities and the calculations from the sta- 
tistical multifragmentation model (lines). The regions 
enclosed between the solid lines correspond to the SMM 
calculation obtained from the secondary fragment yield 
distribution, and those between the dotted lines corre- 
spond to the densities obtained from the primary frag- 
ment yield distribution. The width of the two regions 
represents an upper limit to the sensitivity of the cal- 
culation to the assumed source size. The error bars for 
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the excitation energies corresponds to two different equa- 
tions of state used in the BNV calculation, namely, the 
asy-stiff and the asy-soft equation of state. It is observed 
that the statistical multifragmentation model calculation 
from the secondary fragment yield distribution explains 
the observed data quite well. The calculated densities 
from the primary and the secondary fragment distribu- 
tion are not very different for the 58 Fe + 58 Ni reaction. 
But those for the 58 Fe + 58 Fe reaction, which has a higher 
neutron-to-proton (N/Z) content, are significantly differ- 
ent. The difference is small at lower excitation energy and 
becomes larger with increasing excitation energy. The 
comparison shows that the primary fragment yield dis- 
tribution for systems with large neutron-to-proton (N/Z) 
ratio and initial excitation energy are strongly affected 
by the secondary de-excitation process. With increas- 
ing excitation energy and neutron-to-proton ratio of the 
fragmenting system, the primary fragments are highly 
neutron rich and produced with large internal excitation 
energy which quickly undergoes secondary de-excitation 
leading to a decrease in the width of the yield distribution 
as observed in Figure 1. Such an effect has also been ob- 
served for fragments produced in a dynamical stochastic 
mean field calculation 0, [2(j • 



V. THE REDUCED NEUTRON DENSITY AND 
THE SYMMETRY ENERGY 

In the following, we extract the symmetry energy of 
the primary fragments from the reduced neutron density 
and study their dependence on the excitation energy of 
the fragmenting system. From Eq. 1 the ratio of the 
fragment yield from two different reactions can be written 



Y(N, Z) 
Y m + Ni (N, Z) 



C 



Pn 



Ni+Ni 
Pn 



N 



Pp 



Ni+Ni 
Pp 



= cf n p z v 

(4) 

Where C is the overall normalization factor and p n , 
p p are the relative reduced neutron and proton densities. 
The relative reduced neutron and proton densities are 
related to the isoscaling parameters a = Ap n /T and j3 
= Afip/T through a relationship given by e a = p n and e@ 
— pp. The quantities Ap n and Ap: p are the difference in 
the neutron and the proton chemical potential between 
the two reactions. The parameter a is related to the 
symmetry energy part of the fragment binding energy 
C sym , and given as (see |(|), 



aT 



IP |iiL_£l 
I A\ A\ 



(5) 



where Z\, Ai and Z 2 , A 2 are the charge and the 
mass numbers of the fragmenting systems. Figure 5 (top 
panel) shows the experimentally (symbols) determined 
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FIG. 5: Isoscaling parameter (top panel) and symmetry en- 
ergy (bottom panel) as a function of excitation energy of the 
fragmenting source for the 58 Fe + 58 Ni and 58 Fe + 58 Fe re- 
actions. (Top panel) The dashed and the solid curves are the 
calculated values of a from the primary and secondary frag- 
ment yield distributions using SMM, respectively. (Bottom 
panel) Symmetry energy calculated using the primary frag- 
ment yield distribution. 



isoscaling parameter a, from the reduced neutron den- 
sity as a function of excitation energy of the fragmenting 
system. Also shown in the figure are the calculated val- 
ues of the scaling parameter from the primary (dashed 
curve) and the secondary (solid curve) fragment yield dis- 
tribution using the statistical multifragmentation model. 
Within the uncertainties in the input parameters of the 
SMM calculation, the decrease in the a values with in- 
creasing excitation energy is well explained by the calcu- 
lation once the secondary de-excitation effect in the pri- 
mary distribution is taken into account. Also, the calcu- 
lated values from the secondary fragment distribution are 
not very different from those calculated from the primary 
fragment distributions for the lower excitation energies. 
However, with increasing excitation energy and isospin 
(N/Z) of the system they differ significantly. More specif- 
ically, the scaling parameter a for the neutron-rich pri- 
mary fragments produced at excitation energies greater 
than 5-6 MeV/nucleon are strongly affected by the sec- 
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ondary decay effect. 

Having demonstrated the fact that the observed scal- 
ing parameter is well reproduced by the statistical multi- 
fragmentation model and that a high degree of thermal- 
ization is expected during the fragment formation, we 
have attempted to obtain an estimate of the symmetry 
energy C sym of the primary fragments using equation 5. 
The temperature for the estimate was determined using 
the double isotope ratio of the fragment yield Q . Since 
the determination of the temperature from the fragment 
yield is also affected by the secondary de-excitation, a 
maximum correction of 70% was applied to be consis- 
tent with the previous determination of temperature for 
similar reactions (see [2l| ). This correction results in a 
temperature value of 6 - 7 MeV over the energy range 
studied. The difference {Z/A)\ - {ZfA)\ of the frag- 
menting system was estimated at t = 50 fm/c using the 
dynamical model calculations [ToL fill and was about 
a 3% lower than the difference in the initial Z/A of the 
systems. Within these uncertainties, one obtains C aym 
between 15 - 20 MeV over the excitation energy range 
studied. The estimated values from both the systems, 
58 Fe + 58 Fe and 58 Fe + 58 Ni, are comparable to each 
other and shown by the thin and thick dashed lines in 
Fig. 5 (bottom panel) . These observed values are signif- 
icantly lower than the standard value of 25 MeV (dot- 
ted line) often assumed for a stable and unexcited nu- 
clei. Furthermore, the symmetry energy is also observed 
to decrease slowly with increasing excitation energy. It 
should be mentioned that the estimated value of the sym- 
metry energy is sensitive to the corrections assumed for 
the {Z/A)\ - (Z/A)l in Eq. 5. A lower value of the 
{Z/A)\ — (ZjA)\ could raise the value of symmetry en- 
ergy. However, from the present observation it would 
require a correction of more than 70% at the highest ex- 
citation energy to reproduce the standard value of 25 
MeV. Such an increase in not supported by any dynam- 
ical calculations and appears unlikely. It thus appears 
impossible to reproduce the standard value of the sym- 
metry energy by any means. The data along with the 
calculations thereby suggests that the symmetry energy 
is below 20 MeV and well below the standard value. The 
reduced value of the symmetry energy shows that the 
primary fragments are not only excited and neutron rich 
but also expand to a reduced density when formed. For 
higher excitation energies the fragments appear to ex- 
pand to even lower symmetry energies. A self consistent 
check of the SMM calculation shows strong variations 



in the properties of the primary fragments when C sym 
is varied in the range 15 - 25 MeV. Recently @, it has 
been shown that neutron rich hot nuclei can be produced 
in stellar matter between the protoneutron star and the 
shock front in type II supernova explosion. In particular, 
it has been shown that a slight decrease in the symme- 
try energy coefficient can shift the mass distributions to 
higher masses. This property of hot nuclei could thus be 
interesting to investigate for understanding the relative 
abundance of elements in the core collapse supernovae 
explosion. 



VI. CONCLUSION 

In conclusion, we have studied the isotopic properties 
of the primary and secondary fragment yield distribution 
using the reduced nucleon densities in multfragmentation 
of 58 Ni + 58 Ni, 58 Fe + 58 Ni and 58 Fe + 58 Fe reactions at 
30, 40 and 47 MeV/nucleon. It is shown that the primary 
fragment yield distribution for system with large excita- 
tion energy and neutron-to-proton ratio are significantly 
affected by the secondary de-excitation effects. The effect 
is small at lower excitation energy and smaller neutron- 
to-proton ratio and becomes large with increasing excita- 
tion and neutron-to-proton ratio. The symmetry energy 
for the fragments deduced from the scaling parameter 
after correcting for the secondary decay effect is consid- 
erably lower than the standard value at saturation den- 
sity. Furthermore, the symmetry energy decreases with 
increasing excitation energy indicating that the fragment 
nuclei are hot and expand to a reduced density when 
formed. The present observation could be important in 
the understanding of the mechanism of the core collapse 
supernovae explosions and the relative abundance of the 
produced elements. 
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